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ABSTRACT 

High spatial and sensitivity images of the Luminous Blue Variable 
IRAS 18576+0341 were obtained using the mid infrared imager VISIR at the 
Very Large Telescope and the Very Large Array interferometer. The resulting 
mid-infrared continuum maps show a similar clumpy and approximately circular 
symmetric nebula, which contrasts sharply with the asymmetry that characterizes 
the ionized component of the envelope, as evidenced from the radio and [Ne n] line 
images obtained with comparable spatial resolution. In particular, there is excellent 
overall agreement between the 12.8 //m map and the radio images, consistent with 
free-free emission from circumstellar ionized material surrounding a central stellar 
wind. 

The color temperature and optical depth maps obtained from mid-infrared images 
show only slight fluctuations, suggesting quite uniform dust characteristics over the 
dust shell. 

We explore various possibilities to understand the cause of the diff'erent morphol- 
ogy of the dusty and gaseous component of the circumstellar envelope which are com- 
patible with the observations. 

Subject headings: circumstellar matter — infrared: stars — stars: early-type -stars: individual 
(IRAS 18576+0341) —stars: winds, outflows 
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1. Introduction 



During the latest phases of their evolution, high mass stars experience significant mass loss, 
both through strong stellar winds and eruptive events, that drives these objects out of the main 
sequence through short and poorly constrained instability phases. The class of Luminous Blue 
Variables (LBVs) consists of luminous (log(L/Lo > 5.5) and massive stars that are believed to 
go through a short but violent t ransition from the main sequence towards the Wolf-Rayet stage 



(|Humphreys & Davidson 



1994|) . although thei r link to other advanced evolutionary ph ases of 



massive stars such as supemovae is still open (IBarlow et al. 



2005 



Kotak & Vink 



20061). During 



their evolution, LBVs lose a huge quantity of mass from their original envelope, leading to the 
formation of extended circumstellar nebulae (LBVN), whose masses are thought to be on the 
order of a few solar masses. In the galactic context, despite the relative small number of confirmed 
and candidate members, LBVs represent a significant stage of the cosmic cycle of matter, as 
they contribute to the ch emical enrichment of interstellar medium in dust and heavy ele ments 



( Smith & Qwocki 



2006|) . The census of Galactic LBVs performed by Clark et al. 



(l2005h reported 



12 effective members and 23 candidates, but new candidates have recently been identified via 
the mid-infrared i mages of the Galac tic plane obtained with the Multiband I maging Photomete r 



for Spitzer (MIPS; 



Rieke et al. 



2004 and the Infrared Array Camera (IRAC: 



onboard the Spitzer Space Telescope dWachter et al 



20101 : 



Gvaramadze et al 



Fazio et al 



2010a 



2004) 



The overall scenario of LBV evolution is currently under debate and the nebula formation 
mechanism, i. e. intensity, duration and geometry of the mass-loss events, is not well established. 
The study of the gas and dust content of the nebula which surrounds the star may provide 
significant clues that will enable us to understand the mass-loss processes that govern its evolution, 
as well as the drast i c changes that take place in the immediate circumstellar environment 



dUmana et al 



20091 : 



Jimenez et al. 



2010|) . By exploring the characteristic of the mid infrared 



continuum emission, which traces the warm dust, and of the recombination lines and radio 
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continuum emission, which trace the ionized gas, we can study the relationship between these 
different components which coexist in the stellar ejecta with the aim of gaining valuable insights 
about the formation and shaping mechanisms at work in the circumstellar envelope. In particular, 
radio observations offer the additional possibility to reveal the ionized gas well inside the dusty 
envelope and, therefore, allow to detennine in grea t detail its spatial distribution without suffering 



from intrinsic extinction (Umana et al. 



2005L 



2010 ). 



The infrared source IRAS 18 576+0341 was originally classified as a Planetary Nebula 



on the basis of its in frared colors (|Garcia-Lario et al, 



(Condon et al. 



object (|Ueta et al 



1997h and of the observed 1.4 GHz flux 



1999h. Su bsequently, the detection of mid-infrared extended emission around this 



infrared variability (lUeta et al. 



200 ih. together with observations of strong photometric and spectros copic 



2001 



Pasquali & Comeron 



20021: 



Clark et al. 



2001 



200% and its 



position in the HR diagram led to a definitive classification of IRAS 18576+0341 as a LBV. 



On the basis of radiative transfer calculations. 



Ueta et al, 



(|200l[) derived a log = 6.4 and a 



temperature T^ff = 15 + 9 kK for the central star, which appears s urrounded by a neb ula of gas and 



Clark et al. 



(2002), in a detailed 



dust characterized by a possible C/0 mixed chemistry. Recently, 
analysis of long term monitoring of the near infrared photometric and spectroscopic properties of 
IRAS 18576+0341, showed that it has been highly variable over the last 20 years. Near-infrared 
photometric variability greater than one magnitude has been observed with moderate variation in 
temperature (AT^, < 4.5kK), which is consistent with variations of stellar radius. 

Multi-frequency VLA observations have revealed an extended, asymmetric, and quite 
structured nebula and revealed the central core of the LBV, w hose radio propertie s are typical of 



a stellar wind, with a current mass-l oss of 3.7 x 
consistent with those determined by IClark et al. 



OJMq yr ^ (lUmana et al. 



20051). This value is 



(|2009r) from the NLTE atmospheric code, which 
indicates a variation of the mass loss rate from 1.2 x 10""^ to 5.2 x lO^^M© yr'^ between 2002 and 



2006, associated with a decrement of Teff from 15.5 to 1 1 .0 kK 
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The present work is a part of an effort to better understand the physical conditions in a 
sample of Galactic LBV nebulae, conducted by using state of the art instruments available at mid 
infrared and radio wavelengths. We present the results obtained for the LBV IRAS 18576+0341 
by combining mid infrared and radio maps having comparable spatial resolution, obtained using 
VISIR at Very Large Telescope (VLT) and Very Large Array (VLA). The images allow us to study 
the detailed morphology of the circumstellar envelope, in order to un derstand the origin of the 



apparently contradictory spatial structures seen in the resulting maps (|Umana et al. 



2010 ) 



2. Observations and Data Reduction 

2.1. VISIR/VLT Observations 

To map the dust distribution in the LBVN associated with IRAS 18576+0341, we have 
obtained high angular resolution and high sensitivity N and Q-band images with VISIR 
([Lagage et al. 



2004|) . the VLT imager for the mid-IR mounted at the Cassegrain focus of the VLT 
Unit 3 telescope (MELffAN). The observations were carried out on 2007 July 19 and 26 through 
the PAH2 filters, centered on known R\H features {Ac = 1 1.26/im, AA = 0.59/zm) and adjacent 
continuum (R\H22, Ac = 11.88/im, AA = 0.37/zm), and in Ql {Ac = 17.65/im, AA = 0.83//m). 
We also used the Nell filter {Ac = 12.80/zm, AA = 0.21/im) and its adjacent continuum (Nelli 
Ac = 13.03 jum, AA = 0.22 fim) as a tracer of ionized gas to compare with the radio images. All the 
observations were performed under very good and stable weather, with an optical seeing of about 
0.7". The target was observed at airmasses ranging between 1.2 and 1.4. We used a fixed pixel 
scale of 0.075", resulting in a FOV of 19.2 x 19.2". The standard chopping/nodding technique 
was adopted for subtraction of the sky background as well as the telescope's thermal emission: 
secondary mirror chopping was performed in the North-South direction, with a chop throw of 9" 
and telescope nodding was applied in the opposite direction with equal amplitude. To further 
improve the image quality, a random jitter pattern with a maximum throw of 3" was superimposed 
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on the nodding sequence. A complete log of observations is summarized in Table 1 . Data were 
reduced following the standard VISIR pipeline (version 1.7.0), consisting of co-adding of frames 
after flat fielding and removal of bad pixels. The chopped and nodded images were then combined 
to make one image in each filter. Three mid-IR standards were observed just before and after the 
target acquisition to flux calibrate the data. Such observations were also used to determine the 
FWHM of standard star images and thus to derive the actual angular resolution of our final maps. 
Values of the FWHM for each filter are reported in Table 1 . 



2.2. VLA observations 



We observed IRAS 18576+0341 on the 2004 October 12 at 4.8 GHz (6cm) and 1.4 GHz (20 
cm), with a bandwidth of 100 MHz using the VLA in A configuration, providing a typical beam 
size of ~0.36" and xl3" respectively. For all the observed bands, 1824+107 was used as phase 
calibrator and the flux density scale was determined by observing 3C 48. 

The data processing was performed using the standard programs of the NRAO Astronomical 



Image Processing System (AIPS). We compared the images made from t 



le A-c onfiguration data 



2005 ) to confirm that 



with our earlier images obtained using the C configuration dUmana et al. 
long term variability was not present in our data. 

After the calibration process, the A-configuration wv-datasets were combined with the earlier 
o bservations pe r forme d at the same frequencies using the VLA in C configuration and presented 



in 



Umana et al. 



(I2005|) and images have been produced at both wavelengths. This allowed us 
to obtain high resolution maps of the entire nebula and, in the same time, fully recover all the 
extended emission. The mapping process was performed by using AIPS task IMAGR with a 
variety of weighting schemes and the dirty maps were CLEANed down as close as possible to the 
theoretical noise. The final uniform-weighted (ROBUST -5) maps have a synthesized beam of 
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a'38xa'36 and l'.'39xl'.'27, and rms noise of 0.04 mJy beam"' and 0.1 mJy heam'^ at 6 cm and 
20 cm, respectively. The noise level (rms) in the maps was estimated by analyzing an area on the 
map (using task IMEAN), whose dimension is of the order of more than 100 O^y^ away from the 
phase center and free from evident radio sources. 



3. Results 



3.1. The dust shell 



All five VISIR images, shown in Figure [H have a similar morphology, with an extended 
circumstellar envelope surrounding a central source that is clearly visible until almost 13 fim. In 
Figure [U we show the images obtained in the filters PAH22, Nen2 and Ql. From all the filters we 
derived a shell diameter o f approximately 7", corresponding to a size of 0.35 pc at a distance of 



10 kpc (IClark et al 



2009h 



We performed aperture photometry of the central object in each image by using the IDL 
procedure ATV, while the nebula contribution at each wavelength has been determined by 
integrating the flux emitted in a selected area of 9" x 9". The photonietric results are summarized 



Hrivnak et al. 



(|2000|) pointed out a strong 



in Table [31 On the basis of their ISO SWS spectrum, 
difference between the IRAS and the simulated ISO photometry, concluding that the dusty 
envelope around IRAS 18576+0341 should have been very extended with some of the flux 
observed by IRAS falling outside of the largest ISO aperture (> 30"). However, we fully recover 
the IRAS 12 /um flux density (F60;jm=58.48 mJy), implying that, at least at 12 yum, all of the 
emission is from the compact 7" source. 



200 ih that discovered the 



With respect to the previous mid-infrared images (|Ueta et al.l 
presence of a roughly circularly symmetric nebula around IRAS 18576+0341, our higher 
resolution and sensitivity allow us to probe finer details of the dust distribution. The brightness 
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of the shell varies around its circumference, with the major contrib ution to the flux c oming from 



Uetaetal 



( 2001 ) as an arc 



the north-west part of the nebula. The condensations, detected by 
located in the north side of the shell, now shows up as a clumpy structure with two major peaks 
located at 50 degrees in the East-North direction and abou t 1" from the cen tral star. It appears that 



UetaetaL 



(|200lh were an effect of the 



the northern and the southern emission peaks reported by 
limited sensitivity and spatial resolution of their maps. From our VISIR maps, no hints for the 
claimed optically thin edge-on dust torus surrounding the central object are evident, but only a 
suggestion of a possible spherical inner shell whose edges are strongly structured. 

It is possible to isolate the PAH feature at 1 1 .26 yum by subtracting the image obtained at the 
adjacent continuum (PAH22) from the the PAH2 image. We have performed this subtraction after 
registering the images on the c entral object and normalizing them to the continuum as derived 



from a fit to the ISO spectrum dSloan et al 



20021). The difference image is shown in Figure [21 The 



1 1 .26 yum PAH emission appears to be located in the inner shell with an higher concentration in 
the north-west part of the ring where the most of the thermal dust continuum emission originates. 

The [Ne II] emission, which should trace the ionized gas, can be isolated following the 
same procedure, i. e., by subtracting the adjacent continuum image at 13.03 /zm (Nell2) from the 
12.80 /zm image. The striking result, shown in Figure [3l demonstrates that the [Ne II] emission 
has a very different distribution compared to the dust and PAH tracers. As [Ne II] emission is 
a ionized gas tracer this result suggests that significant difference exist between the dusty and 
ionized stellar ejecta (sec. 13.21) . 

Assuming that the extinction is constant across the images and that the flux that we have 
mapped is due to optically thin thermal emission from dust grains, we may use the NeII_2 and Ql 
images to produce a ratio map. In this case the observed intensity is given by 



(1) 



and using the Wien approximation {hv » fcr) for the blackbody expression By(T), the ratio of the 
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intensities at two different frequencies vi and V2 is given by 



The optical depth t,i depends on the chemical composition and on the size a of the grains. 
There is an open co ntroversy on the nature of dust grains present in the nebula surrounding 



(2) 



IRAS 18576+0341 (UetaetaL 



2001 



Hrivnak et al. 



20001) . indicating that it is possibly a 



mixed-chemistry object where feature s due to silicates and those related to carbon (PAH) coexist. 



In this paper, we follow 



Umana et al. 



torn and adopt for = 11.88 and A2 = 17.65 yum, 



taJtaj = 1-05 in the case of graphite and 1.27 in the case of astronomical silicates. 



For each pixel of the map, we derive the dust temperature inverting equation |2] 

T « 1.44x10^- 



In 



■'"-ii V -^2 / 



(3) 



The dust temperature has been computed only at those pixels whose brightness is greater 
than 0.06 Jy arcsec"^, corresponding to 3 cr in the maps. In this derivation, we assume a constant 
temperature along the line of sight for each pixel. 

The resulting dust temperature map is shown in the left panel of Figure HI The temperature 
map can be used to examine dust temperature gradients across the nebula. We see that the color 
temperature values are confined to a relatively narrow range, from ~ 130 to 160 K in case of 
graphite grains. Slightly lower values (~ 10%) are obtained if instead a mixture of silicate grains 
is assumed. 

The higher temperatures are reached at a distance of about 2" from the center. To more 
clearly see the temperature distribution as a function of radius, in the right panel of Figured we 
show radial cuts of the temperature map along diameters obtained at steps of 10°. 

To evaluate the mass of the dust in the envelope, an optical depth map at A = 17.65 fim is 
required. This can be obtained by inverting EqlH using the value of the temperature found in each 
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point of the map and the observed intensity. If the medium is homogeneous, the optical depth is 
given by T = ^dustPdust U where p^ust and / are the the density of the medium and the thickness of 
the envelope, and ^dust is the absorption coefficient per mass unit. 

The mass behind each optical depth map pixel, whose area is (y.'OVS x (y.'OVS , is given by 
Am = Tik X AA, and AA the are a source correspon ding to one map pixel, which is 1.25 x 10^^ cm^ 



assuming a distance of 10 kpc (|Clark et al 



20091). A value of /tdust = 253 and 1082 cm^ g"^ 



at A = 17 fxm has been adopted in case of graphite and silicates, respectively (lUmana et al 



2010h . Integrating over all the pixels with brightness higher than 3 cr, a total dust mass of 
8.6 X 10^ |g (4.3 X lO'^Mp,) is derived in case of graphite, and 1.3 x lO^^g (6.5 x IO^^Mq) in case of 
silicates. Ueta et all ( 2001 ) provide an estimate of total dust mass of ~ O.IMq assuming a toroidal 
morphology for the entire dust emitting region and a radius of ~ 14 '. This value is compatible 
with our results if scaled down to the extent of the mid-infrared emitting dust shell that is resolved 
in our images. Because the mid-IR observations we present only trace the cooler dust component 
in the nebula, our dust mass estimates should be considered as lower limits to the total dust mass 
in the LBVN. 

The dust column density along the line of sight Pdust I (Figure |5l left panel) can be used to 
infer the density pdust inside the nebula. Although its distribution is clumpy, as already noted, with 
a higher concentration of dust in the north-east part, a roughly spherical symmetry around the 
central star can be assumed. Cuts of the column density map along diameters, obtained at steps 
of 10° averaged together, are shown as a thin line in the right panel of Figure [51 We modeled 
the density distribution assuming a radial dependence. For this purpose, we built a 3D matrix, 
sampled with a step corresponding to the pixel size of the VISIR maps (1.12 x 10'^ cm) in all the 
three spatial dimensions. We found that the density must increase as a function of radius from 
the center, although it is not possible to discriminate among several radial dependencies. As an 
example, the radial column density profile obtained assuming pdust °^ r^, with a smooth at the 
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edge, is shown as a thick line in Figure [51 right panel. 

This analysis shows that in the inner part of the nebula the density is low and it increases 
strongly with the distance from the star. This could indicate: 

a) an episode of strong mass loss at a previous time; 

b) the dust is processed by the UV radiation field of the star in the inner part of the nebula. 

c) a combination of the two. 

The maximum of pdust is at 3" from the center, corresponding to a linear distance of 0.15 pc; 



assum ing the same expansion velocity Vg ~ 70 km s measured from the Bry line ((Clark et al 



20091), the mass loss episode should have occurred about 2 000 years ago. 



3.2. Ionized gas 

The continuum emission associated with IRAS 18576+0341 at 5 GHz and 1.4 GHz is shown 
in Figure [6l The spatial resolution achieved in the final image is about V.'3 in the L Band and 0'.'4 
in the C Band. The images confirm the structure of the source as detected at higher frequencies 



( Umana et al. 



2005|) : a compact, slightly resolved core, most likely associated with the stellar 
wind from the central object, surrounded by an extended ionized nebula with a remarkably 
similar shape and extension at all wavelengths. As expected, the combined eff"ect of reduced 
spatial resolution and of intrinsic contribution at low frequency prevent us from detecting the core 
component and discerning it from the extended nebula in the 20 cm map; on the contrary, the 
higher spatial resolution in the 6 cm map allows us to clearly resolve the compact source. 

The total radio flux densities associated with the entire source have been derived by summing 
the contribution of all the pixels contained in the radio source (AIPS task IMSTAT), while the 
flux density of the central object object has been obtained by fitting a two dimensional Gaussian 
brightness distribution to the map (AIPS task JMFIT). The results are reported in Tabled 
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To better separate the core component from the remaining diffuse emission, the flux density 
from the compact core at 6 cm has been determined by fitting a two dimensional Gaussian 
brightness distribution to the map obtained with the A configuration data alone. The resulting 
flux density of 1.3+0.1 is in accordance with the flux density extra polated from measur ements at 



higher frequencies assuming a stellar wind spectrum, as reported in 



Umana et al 



mm. 



The final images in Figure [6] show that the radio emission from the circumstellar ionized 
nebula is asymmetric with respect the central core, extending out in the northwest direction, with 
the brightest radio region located about 3'.' 6 to the north of the central core in the 6 cm map. The 
diffuse emission is remarkably similar to that detected at all the other observing frequencies, with 
similar extent of about 7". The southern component, which contributes about 25% of the total 
mid- infrared emission, is not present in any of the radio maps. A ce rtain degree o f asymmetry is 



also present in the radio emission of other LBVNs, such as HR Car (jWhite 



(|Lang et al 



2005 



20001) and Pistol Star 



19991). This may indicate that extremely asymmetric mass loss can be quite 



common in LB Vs. 

The overall morphology of the radio nebula agrees very closely with the spatial distribution 
of the [Ne II] line emission. The contours of the [Ne II] map, superimposed on a gray-scale 
image of the 6 cm data, is shown in Figure [3l It is evident how the emission in the [Ne II] line 
follows quite faithfully that of the radio continuum. Such resemblance between the two maps is 
unsurprising, as both are tracers of the thermal ionized gas. 

To search for the existence of denser regions with different opacity through the nebula, we 
computed the spectral index map by comparing the 6 cm with the 0.7 cm image obtained in the 
first epoch. The images have comparable spatial resolution. We computed the spectral index, 
defined by the relation S(v) « v", using the maps produced with the same beam and cellsize for 
each pixel with a flux density of at least three times the rms noise at both the wavelengths. The 
resulting map, shown in Figure |7l shows evidence of the stellar wind component, characterized by 
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spectral index of a ~ 0.8, and the extended halo having a much flatter radio spectrum (a ~ -0.1), 
which is typical of free-free thermal emission. However, the image shows some patchiness 
probably associated with local density clumps. We note the presence of a second component with 
spectral index of a « 0.9 located about I'.'S to the south of the central core, in the direction of the 
southern infrared emission peak in the dust nebula. Such a component, which is clearly visible in 



the radio map reported in Figure lb in 



Umana et al. 



(|2005h . is probably due to the presence of a 



dust condensation whose thermal emission contributes to the flux observed at 7 mm. 



4. Discussion 

The most striking result of our observations is the clearly dissimilar morphology of dusty and 
gaseous components of the nebula, which seem to have quite different spatial distribution. This is 
immediately apparent in Figure [81 which shows the composite image obtained by superposition 
of the 6 cm and the 17 fim images. The ionized nebula, in fact, extends beyond the dust shell in 
the N-E direction, and the radio and [Ne II] maps do not contain any indication of the symmetric 
shape seen in the mid-infrared continuum maps. 

Diff"erent scenarios can be envisioned that could account for all the discussed findings and 
provide a complete picture of the source. As the continuum infrared maps seem to exclude 
an asymmetry in the mass loss, the asymmetry in the distribution of the ionized gas could be 
ascribed to one of the following possibilities: (a) direct photo-ionization from an external source, 
(b) variation in the density of the interstellar medium into which the ejected material is expanding, 
or (c) UV radiation from the central star leaking through a hole in the opacity of the circumstellar 
nebula. 

The hypothesis of direct photo- ionization canno t be ruled out because of the high visual 



extinction of roughly Av/ ~ 28 mag (|Ueta et al. 



200 ih towards IRAS 18576+0341. The surveys 
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available at this moment, in fact, do not have the sensitivity needed to detect a main sequence star 
of spectral types O to B as it would appear at a distance of 10 kpc in this direction. For the same 
reason, not even the inspection of our radio maps allow us to exclude the presence of such kind 
of source, although no radio emission has been detected in our map around the location of the 
nebula. On the basis of the analysis from Panagia & Felli (1975), in fact, the flux density expected 
at 6 cm from a spherically symmetric, isothermal and fully ionized stellar wind at a distance of 10 
kpc is about 1.6x10^ mJy, assuming a stellar temperature of 10"^ K, a terminal wind velocity of 
10^ km s"' and a mass loss rate of 10"^ Mq yr~^. The resulting flux is thus very close to the noise 
of the radio map. 

On the other hand, based on at the near- infrared images from the Two Micron All Sky Survey 
(2MASS) toward the position of IRAS 18576+0341, we found that north-east rim of the ionized 
cloud is facing a source designed a s 2MASS 19001104+0345511 in the 2MASS All-Sky Catalog 



of Point Sources (Cutrietal 



120030 and located at position a = 19''00'"ir.05 6 = +03°45'51".2. 
This is a distance of 4'.'7 from the central star of IRAS 18576+0341, just in the direction where 
the nebula ionization appears to be strongest. However, based on the above considerations and the 
near-IR colours derived from the de-reddened 2MASS magnitudes, we exclude the possibility that 
such a source could be an early type main sequence or supergiant star able to provide the energy 
needed for the nebula ionization. 

We can also consider the hypothesis that the ionized gas surrounding the LBV originates from 
the interaction of the mass outflow from the central star with a high density medium surrounding 
the star, such as clumps or perhaps circumstellar material. We can make an o rder-of-magnitude 



Curiel et al. 



(Il989h . 



check of the plausibility of this scenario on the basis of the theory outlined by 
who modeled the shock-ionization scenario and derived the radio continuum emission under 
optically thin conditions. In their formulation, the effective optical depth of the emitting region is 
directly connected to the outflow characteristics, and a simple relation exists between the radio 
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flux density in the optically thin regime and the wind momentum rate: 
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In this formula, 5 y is the radio flux observed at frequency v, rj represents the fraction of the 
stellar wind that is shocked and produces the observed radio continuum emission, M is the mass 
loss rate, Voo is the terminal velocity of the wind, T is the wind temperature and d the distance 
to the source. Adopting rj x 0.3 and a temperature of lO'^K, the density flux of about 100 mJy 
observed at 6 cm can be provided by isotropic mass loss rate ranging from 7x10"^ to 3 xlO^ 
Moyr~' for an assumed wind velocity varying from 500 to 200 kms"\ which are parameters 



values co mparab 



eruptions (ISmith 



e with those derived for IRAS 18676+0341 and other LBVs observed giant 



2005 



Smith & Hartigan 



20061: 



Clark et al. 



20091). Due to the many unknown 



parameters this estimate is only suggestive, but it shows that with typical parameters the shock 
emission model can produce the observed flux. 

On the other hand, in absence of an external factor (source), the ionization of the circumstellar 
nebula could be explained exclusively by the direct photo-ionization from the central object. In 
this case, as derived in Paper I, the Lyman continuum flux of a BO-BO. 5 supergiant can account 
for the measured radio flux. This scenario, however, involves the presence of some kind of 
inhomogeneity (asymmetry, anisotropics) in the circumstellar material that has not been observed 
in our mid-IR images. It is, in fact, difficult to understand the departures from symmetry of the 
radio and [Ne II] images except in terms of density or optical depth variations in the cloud. 

The very different spatial distribution of circumstellar dust and ionized gas might reflect an 
anisotropy in the same mass loss event that could have excavated the mid-IR emitting dust, or may 
indicate that distinct ejection events have been coUimated differently by the interaction with the 
circumstellar material ejected in previous phases. Such interactions could have affected the grain 
size distribution throughout the nebula, causing an anisotropy in the opacity of the dust grains 
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which are competing with the gas in the absorption the UV radiation (|van HofF et al 



20041) . 



5. Conclusions 

We have presented new high-resolution mid-infrared and radio imaging observations of 
IRAS 18576-1-0341. These observations provide an unprecedented view of the detailed spatial 
structure of the circumstellar envelope surrounding the central massive star, and highlight the 
different distribution of dust and gas component in the nebula. These data allow us to perform a 
spatially resolved analysis of dust temperature and optical depth, and to derive the high resolution 
spectral index map of the observed radio emission. 

We found that the dust properties seem to be essentially the same everywhere in the nebula 
and the spectral index map is consistent with optically thin free-free radiation, except for the 
wind component from the central star. We thus conclude that the dumpiness observed at all 
wavelengths has to be ascribed to local density enhancements. 

We also described the possible scenarios that could account for the asymmetry in the spatial 
distribution of the ionized material. We suggest that the gas ionization could be caused by 
UV-radiation from an external source or by the impinging of the mass outflow upon an obstacle 
like a denser circumstellar medium. Alternatively, the asymmetry could be ascribed to UV 
photons leaking out from the central object thought a hole in the circumstellar envelope. 

Although we are unable to rule out any of the discussed frameworks, it is worthwhile to note 
that the last scenario conflicts with the relative minor spatial variations of dust temperature and 
optical depth resulting from the analysis of our mid-IR images. 
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Fig. 1 . — VISIR images of IRAS 18576+0341 and its surrounding nebula taken through the filters 
PAH22,NeIl2andQl. 




Fig. 2. — Continuum subtracted 1 1.26 yum image of IRAS 18576+0341, showing the R^H feature 
emission. 
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Fig. 3. — Continuum subtracted 12.8 /im contour map of IRAS 18576+0341, showing the [Ne 11] 
line emission, superimposed to the 6 cm image. 

Table 1 : Observations Log 



Filter 


Date 


UT 


Integration time 


Airmass 


FWHM 








(s) 




// 


PAH2 


2007-07-26 


04:56 


700 


1.2 


0.27 


PAH22 


2007-07-26 


05:15 


700 


1.3 


0.28 


Nen 


2007-07-19 


06:21 


700 


1.4 


0.30 


Nell2 


2007-07-19 


06:01 


700 


1.4 


0.32 


Ql 


2007-07-26 


05:48 


180 


1.4 


0.38 
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Fig. 4. — Left panel: map of dust color temperature as determined from Ai = 11.88 and A2 = 
17.65 //m maps. For the dust a Carbon chemistry has been assumed. Regions of warmer and cooler 
emission are seen as darker and lighter shades respectively. Right: Temperature profiles along 
diameters at step of 10°; here the central part with the star is not considered; the higher values are 
reached at about 2" from the center. 
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Fig. 5. — Left: Map of dust the optical depth. The greatest concentration of mass is located in the 
north-east part of the nebula. Right: Column density profile: the average of the observed diagonal 
cuts is shown as dashed line; here the central part with the star is not considered; a profile obtained 
with the density increasing outward is shown as a continuous line. 
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Fig. 6. — Uniform weighted radio maps of IRAS 1 8576+034 1 at 6 and 20 cm, obtained by merging 
A and C VLA configuration data. 




Fig. 7. — Spectral index distribution between 6cm and 0.7cm. Contours levels are -0.2 to 0.7 in 
steps of 0.1 
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Fig. 8. — Composite image of IRAS 18576+0341 obtained by superposition of radio images at 
6 cm (red) and 17.65 fxm (blue). See the electronic edition of the Journal for a color version of this 
figure. 

Table 2: Observed Mid Infrared Flux 



Filter 


p 

^ star 


Fneb 


/im 


Jy 


Jy 


11.26 


0.12+0.02 


41.2+0.2 


11.88 


0.15 +0.02 


51.9 +0.2 


12.80 


0.19+0.03 


78.8 +0.1 


13.04 


0.17+0.03 


74.7+0.1 


17.65 




224.7+0.1 
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Table 3: Observed radio flux 



Filter 


Fstar 


Fneb 


cm 


mJy 


mJy 


6 


1.3+0.04 


103+2 


20 




92+2 



